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Abstract 

Models of dark energy are conveniently characterized by the equation-of-state parameter w = p/p, 
where p is the energy density and p is the pressure. Imposing the Dominant Energy Condition, 
which guarantees stability of the theory, implies that w > — 1. Nevertheless, it is conceivable 
that a well-defined model could (perhaps temporarily) have w < — 1 , and indeed such models 
have been proposed. We study the stability of dynamical models exhibiting w < — 1 by virtue of 
a negative kinetic term. Although naively unstable, we explore the possibility that these models 
might be phenomenologically viable if thought of as effective field theories valid only up to a certain 
momentum cutoff. Under our most optimistic assumptions, we argue that the instability timescale 
can be greater than the age of the universe, but only if the cutoff is at or below 100 MeV. We 
conclude that it is difficult, although not necessarily impossible, to construct viable models of 
dark energy with w < — 1; observers should keep an open mind, but the burden is on theorists to 
demonstrate that any proposed new models are not ruled out by rapid vacuum decay. 
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I. INTRODUCTION 



Cosmological observations strongly indicate that the universe is dominated by a smoothly 
distributed, slowly varying dark energy component ; for reviews see |3j, |4| . ) The simplest 
candidate for such a source is vacuum energy, or the cosmological constant, characterized by 
a pressure equal in magnitude and opposite in sign to the energy density: 

Pvac Pvac • (1) 

While vacuum energy is strictly constant throughout space and time, it is also worthwhile 
to consider dynamical candidates for the dark energy. A convenient parameterization of the 
recent behavior of any such candidate comes from generalizing the vacuum-energy equation 
of state to 

P = wp , (2) 

which should be thought of as a phenomenological relation reflecting the current amount of 
pressure and energy density in the dark energy. In particular, the equation-of-state parameter 
w = p/p is not necessarily constant. However, given that there are an uncountable number 
of conceivable behaviors for the dark energy, a simple relation such as (J2J) is a useful way to 
characterize its current state. 

The equation-of-state parameter is connected directly to the evolution of the energy den- 
sity, and thus to the expansion of the universe. From the conservation-of-energy equation 
for a component pi in a Robertson- Walker cosmology with scale factor a(t) and Hubble 
parameter H = a/ a, 

pi = -3H(pi+pi) , (3) 
it follows that this component evolves with the scale factor as 

— = -3{l + w i ). (4) 

a In a 

We notice in particular that vacuum energy remains constant, while the energy density 
would actually increase as the universe expands if Wi < — 1. The Friedmann equations may 
be written as 

2 8nG k 
H 2 = —P ~ -5 , (5) 
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where k is the spatial curvature, and 
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a 




(6) 



From (jSJ, we see that the universe will accelerate (a > 0) if w < —1/3. (Of course this is 
the effective w of all the energy in the universe; if there is a combination of matter and dark 
energy, the dark energy will have to have a more negative w in order to cause acceleration.) 
From (0) we see that a flat universe dominated by a component with constant w will expand 

as 



unless w = —1, for which the expansion will be exponential. (For w < —1, one should choose 
t < in this expression.) 

What are the possible values w may take? It is hard to make sweeping statements about 
a component of energy about which we know so little. In general relativity, it is conventional 
to restrict the possible energy-momentum tensors by imposing "energy conditions". In j^] 
it was suggested that a reasonable constraint to impose would be the null dominant energy 
condition, or NDEC (see Section |H] for discussion). The physical motivation for a condition 
such as the NDEC is to prevent instability of the vacuum or propagation of energy outside 
the light cone. Applied to an equation of state of the form Q, the NDEC implies w > — 1. 
Thus, pure vacuum energy is a limiting case; any other allowed component would diminish 
in energy as the universe expands. 

Given our ignorance about the nature of the dark energy, it is worth asking whether this 
mysterious substance might actually confound our expectations for a well-behaved energy 
source by violating the NDEC. Given that the dark energy should have positive energy den- 
sity (to account for the necessary density to make the universe flat) and negative pressure 
(to explain the acceleration observed in the supernova data), such a violation would imply 
» < -1. It has bee,, know,, for 8 o,„e time that sueh energy components can occur flQfl. 
Their role as possible dark energy candidates was raised by Caldwell |9(, who referred to 
NDEC-violating sources as "phantom" components, and has been since investigated by sev- 
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(7) 
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eral authors (for some examples see 



limits on w 
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0, ii, is,Q). 



Observational 



21] are conventionally expressed as allowed regions in the w-Qm plane, as- 



suminga flat universe (k = 0, or Om+^x = 1, where X stands for the dark energy). Current 



limits 



2l| . obtained by combining results from cosmic microwave background experiments 



with large scale structure data, the Hubble parameter measurement from the Hubble Space 
Telescope and luminosity measurements of Type la supernovae, give —1.62 < w < —0.74 at 
the 95% confidence level. 

It is straightforward to examine the cosmological consequences of a dark energy compo- 
nent which is strictly constant throughout space, and evolving with any value of w. Any 
physical example of such a component, however, will necessarily have fluctuations (so long 
as w / —1). It is therefore important to determine whether these fluctuations can lead 
to a catastrophic destablization of the vacuum. Unfortunately, because p = wp is a phe- 
nomenological description valid for a certain configuration rather than a true equation of 
state, specifying w is not enough to sensibly discuss the evolution of perturbations, since 
Sp w5p. We must therefore choose a specific model. In particular, the simplest way to 
obtain a phantom component (w < — 1) is to consider a scalar field with negative kinetic 
and gradient energy p, 

1 ._ 1 

(8) 



I(V0) 2 + K(0) 



Fluctuations in this field have a negative energy, and it may be possible for the vacuum to 
decay into a collection of positive-energy and negative-energy particles. If the timescale for 
such an instability is less than the age of the universe, the phantom component would not 
be a viable candidate for dark energy. 

Our goal in this paper is to ask whether phantom components are necessarily plagued 
by vacuum instability, and hence whether observers should take seriously the possibility 
that w < — 1. We will start with a survey of energy conditions and model- independent 
considerations in Section [TTJ In Section HTT1 we describe the cosmology of w < —1 models 
more thoroughly and investigate linear perturbations in a cosmological model, demonstrating 
that it can be compatible with current observations. To investigate stability beyond linear 
order, in Section IIVI we consider a classical toy model of a phantom harmonic oscillator 
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coupled to an ordinary oscillator, and demonstrate numerically that, for small perturbations 
and sufficiently small values of the coupling constant, there exist both stable and unstable 
regions of parameter space. 

The heart of our paper is in Section where we consider the field theory of a phantom 
component coupled to gravity, and calculate the decay rate of a single phantom particle 
into several phantoms plus gravitons. The rate is naively infinite, due to the infinite phase 
space of high- momentum particles. We argue that, considering the phantom Lagrangian 
as an effective theory, the rate may be rendered finite by imposing a momentum cutoff and 
demonstrate that, if we restrict attention to couplings in the potential, then, for a momentum 
cutoff not far below the Planck scale, and for a suitable potential, phantom quanta may be 
stable against decay into gravitons and other particles over a timescale long compared to the 
age of the universe. However, when we include derivative couplings of phantom particles to 
gravitons, we find that such operators can lead to unacceptably short lifetimes for phantom 
particles. 

II. CLASSICAL ENERGY CONDITIONS 

In classical general relativity, without having a specific model for the matter sources, we 
can nevertheless invoke energy conditions which restrict the form of the energy-momentum 
tensor T^ v . In this section we will briefly review these energy conditions, discuss which are 
relevant for cosmology, and compare them to the condition w > — 1. A related discussion 
can be found in |22j ]. 

Each of the energy conditions can be stated in a coordinate-invariant way, in terms of 
T^y and some vector fields of fixed character (timelike/null/spacelike). For purposes of 
physical insight, it is often helpful to consider the case of a perfect fluid, for which the 
energy-momentum tensor takes the form 

T^u = (p + p)U^U» + pg^ , (9) 

where p is the energy density, p the pressure, the fluid four-velocity, and g^ v the metric. 
(More precisely, p and p are the energy density and pressure as measured in the rest frame 
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a) WEC I 




b) NEC 




c) DEC 




d) NDEC I 




e) SEC 



f)w>=-l 





FIG. 1: Shaded regions in the p-p plane are those which obey the designated energy conditions. 
Illustrated are the Weak Energy Condition (WEC), Null Energy Condition (NEC), Dominant 
Energy Condition (DEC), Null Dominant Energy Condition (NDEC), the Strong Energy Condition 
(SEC), and the condition w > —1. Definitions of each condition are found in the text. 



of the fluid, but the shorthand designations are standard.) Our metric signature convention 
throughout this paper is ( — h++). 

The most common energy conditions are the following: 
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• The Weak Energy Condition or WEC states that T^t^f > for all timelike vectors 
t M , or equivalently that p > and p + p > 0. 

• The Null Energy Condition or NEC states that T^i^i" > for all null vectors £ M , or 
equivalently that p + p > 0. 

• The Dominant Energy Condition or DEC includes the WEC (T^tH" > for all 
timelike vectors as well as the additional requirement that T^t^ is a non-spacelike 
vector (i.e., that T^ lv T v \tH x < 0). For a perfect fluid, these conditions together are 
equivalent to the simple requirement that p > \p\. 

• The Null Dominant Energy Condition or NDEC is the DEC for null vectors only: for 
any null vector i^, T^^f > and T^l^ is a non-spacelike vector. The allowed 
density and pressure are the same as for the DEC, except that negative densities are 
allowed so long as p = —p. 

• The Strong Energy Condition or SEC states that T fJiU t^t lJ > ^T x \t a t a for all timelike 
vectors or equivalently that p + p > and p + 3p > 0. 

In Figure ^ we have plotted these conditions as restrictions on allowed regions of the p-p 
plane. We have also plotted the condition w > — 1 for comparison. Note that w > — 1 is not 
equivalent to any of the energy conditions, although it is implied by the WEC, the DEC, 
and the NDEC. 

The different energy conditions are used in different contexts — for example, the WEC and 
SEC are used in singularity theorems — and we will not review them in detail here (see 
for a discussion). Our present concern is to understand under what conditions a hypothetical 

e. For this purpose, the relevant 



23, 



24j . The conservation theorem 



dark energy component would be guaranteed to be stab 
result is the "conservation theorem" of Hawking and Ellis 
invokes the DEC, and uses it to show that energy cannot propagate outside the light cone; 
in particular, if T^ v vanishes on some closed region of a spacelike hypersurface, it will vanish 
everywhere in the future Cauchy development of that region — energy-momentum cannot 
spontaneously appear from nothing. A source obeying the DEC is therefore guaranteed to 
be stable. 
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For cosmological purposes, however, the DEC is somewhat too restrictive, as it excludes 
a negative cosmological constant (which is physically perfectly reasonable, even if it is not 
indicated by the data). This is why [5] advocated use of the NDEC in cosmology, since 
the NDEC is equivalent to the DEC except that negative values of p are allowed so long 
as p = —p. Of course, the less restrictive NDEC invalidates the conservation theorem, as 
a simple example shows. Consider a theory with a negative vacuum energy p vac < 0, and 
an ordinary scalar field if) (not a phantom) with potential V(if)) = ^m?if) 2 . Then both the 
vacuum energy and the scalar field obey the NDEC (although their sum may not, since the 
NDEC is a nonlinear constraint.) Imagine a field configuration with if) = ^J—2p vac /m 2 and 
if> — everywhere throughout space. The energy-momentum tensor of this configuration is 
exactly zero [since V(if>) = — p vac ], but it would instantly begin evolving as if> rolled toward 
the minimum of its potential. There is nothing unphysical about such a situation, however, 
since a stable state is achieved once the field reaches if) = 0. The NDEC, therefore, does not 
guarantee stability with the confidence that the DEC does. However, it seems legitimate to 
ask that dynamical fields obey the DEC, while the cosmological constant is allowed as an 
exception. 

The DEC implies that w > — 1. (Indeed, for cosmological purposes we are interested in 
a source with p > 0; in that case, all of the energy conditions imply w > — 1.) Therefore, 
if we allow for pressures which are less than — p, we cannot guarantee the stability of the 
vacuum. The converse, however, is not true; a phantom component will not necessarily allow 
vacuum decay. In fact, since p = wp is just a convenient parameterization and not strictly- 
speaking an equation of state, the issue of stability is somewhat complicated, as we shall see 
in the next few sections. Nevertheless, it seems very likely that energy sources which violate 
the DEC generally will be unstable, and this is what we will find in the specific example 
considered in Section Our philosophy, therefore, is neither to dismiss the possibility of 
DEC violation on the grounds that we cannot prove stability, nor to blithely accept DEC 
violation on the grounds that we cannot prove instability, but instead to see whether it is 
plausible that the timescale for instability could be sufficiently long so as to be irrelevant for 
practical purposes. 
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III. COSMOLOGICAL EVOLUTION AND PERTURBATIONS 



Consider a flat Robertson- Walker universe with metric 



ds 2 



-dt 2 + a 2 (t)[dx 2 + dy 2 + dz 2 } 



(10) 



for which the Einstein equations are the Friedmann equations (0) and (JBJ). The cosmology 
and fate of a universe containing an energy component with constant w < — 1 are relatively 
simple and have been examined in, for example, JjJ. As an example, consider the case in 
which the universe contains only dust and phantom matter. Then, if the universe ceases to 
be matter-dominated at cosmo logical time t m , then the solution for the scale factor is 



a(t) = a(t m ) 



-w 



(l + w)[- 



2/3(1+™) 



From this expression it is easy to see that phantom matter eventually comes to dominate 
the universe and that, since the Ricci scalar is given by 



R 



4(1 -3w) 



3(1 



w 



t 



w 



1 -2 



(12) 



.1 +Ws 

there is a future curvature singularity at t = wt m / (1 +w). This occurs because, even though 
the energy densities in ordinary types of matter are redshifting away, the energy density in 



phantom matter increases in an expanding universe. Thus, the fate o 



29 



y 



31 



;he universe in these 



32j in w > — 1 dark 



models [27[ may be very different from that expected 
energy models. 

It is, however, simple to construct models of phantom energy in which a future singularity 
is avoided. As we shall see in this section, scalar field models can yield a period of time in 
which the expansion proceeds with w < — 1 and yet settles back to w > — 1 (in our case 
w = — 1) at even later times, thus sidestepping the predictions of w — constant models. 



Consider a scalar field theory with action S = J CJ\g \ d x, and Lagrange density given by 



C^ l -g^{dMdA)-V{4>) ■ 



(13) 



The notable feature of this model is that the sign of the kinetic term is reversed from its 
conventional value [in our conventions the usual expression would be — 1(<90) 2 ]. The equation 
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of motion for becomes 

<j> + 3H<j> - a" 2 V 2 - V\(j>) = , (14) 

where V 2 is the spatial Laplacian, V 2 = <9 2 0+<9 2 0+<9 2 0, and a prime denotes differentiation 
with respect to <fi. Here the energy density p^ and the pressure p^ for a homogeneous <fi field 
are given by 

P<t> = Af + V{(j>) , (15) 

P* = Af - V(<f>) , (16) 
so that the equation-of-state parameter 

p |0 2 + V(0) 
w = - = f-: (17) 

P |0 2 -W) 

satisfies w < — 1. 

We are interested in the cosmological evolution of this model and in the behavior of 
linearized perturbations of the phantom scalar field in the resulting cosmological background. 
As noted in |j , the spectrum of fluctuations of a phantom field evolves similarly to that for 
a quintessence field. Consider metric perturbations in the synchronous gauge, 

ds 2 = -dt 2 + a 2 (c% + h^dx'dx 5 . (18) 

A Fourier mode of the phantom field 

0(t,k) = -= / 0(t,x)e- ik - x d 3 x (19) 

\/27r J 

satisfies the equation of motion 

H k + 3H50 k + (k 2 - V")5<p k = (20) 

where h is the trace of the synchronous gauge metric perturbation hij. The effective mass for 
the perturbation is (k 2 — V") 1 ^ 2 . On large scales one may worry that this effective mass could 
become imaginary for a positive V" . We note, however, that a similar problem exists for a 
canonical scalar field for negative V", and is thus easily avoided by the choice of potential; 
in particular, we may choose a potential with V" < 0. An analysis of perturbations in more 
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FIG. 2: The gaussian potential energy of (|21j) . The phantom scalar will evolve to the top of the 
hill and oscillate around the maximum. 

general non-canonical scalar field models is given in 2^]. In this paper, we will examine by 
explicit calculation the evolution of fluctuations in a specific model. 

Since scalar fields with negative kinetic terms evolve to the maxima of their classical 
potential, we consider a gaussian potential, 

V(4>) = V e-^ 2/a2) , (21) 

where Vb is the overall scale and a is a constant describing the width of the gaussian. The 
potential is represented in Figure 121 We obtained the cosmological evolution by numerically 
solving the equations of motion (jSJ), ©, and (JHJ). The initial conditions were ^initial = Mp 
and initial = 0, where M P = (87TG) -1 / 2 is the reduced Planck mass. Since we require dark 
energy domination at the present epoch in the universe, with Qx ~ 0.7 we choose a = Mp, 
V = 3Mpifg. The results are plotted in figures QUI and 151 

Note that during the initial stages of evolution the field <fi is frozen by the expansion and 
acts as a negligibly small vacuum energy component (with w ~ —1). However, at later times 
the field begins to evolve more rapidly towards the maximum of its potential, the energy 
density in the phantom field becomes cosmo logically dominant, and during this period the 
equation of state parameter is much more negative. Finally, in the very late universe, the 
field comes to rest at the maximum of the potential and a period of w = —1 acceleration 
begins. Since w is no longer less than —1, this ensures that there is no future singularity; 
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FIG. 3: Evolution of the density parameters in radiation (Or), matter (Om)> and the phantom field 
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FIG. 4: Evolution of the phantom field (ft as a function of the scale factor. 
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FIG. 5: Evolution of the equation of state parameter w for the phantom field as a function of the 
scale factor. 

rather, the universe eventually settles into a de Sitter phase. 

Using the formalism for calculating fluctuations in a general matter field developed in 
. one may calculate the effects of the fluctuations in this phantom field on the cosmic 
microwave background radiation (CMB). The resulting power spectrum of CMB fluctuations 
is shown in figure E3 We have not attempted a detailed parameter-fitting to CMB data, but 
it should be clear that our phantom model with the potential (J2*T|) does not predict any 
significant departures from conventional dark-energy scenarios; in particular, there is no 
evidence of dramatic instabilities distorting the power spectrum. (A more detailed study of 
;he effect on the CMB power spectrum of phantom fields with a constant w can be found in 
l^.) However, the formalism for this analysis was based on linear perturbation theory, and 
it is quite plausible that instabilities only become manifest at higher orders. In the following 
sections we address this issue, first through numerical investigation of a model with two 
oscillators, and next by calculating the decay rate of phantom particles into phantoms plus 
gravitons. 
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I 

FIG. 6: Predicted multipole moments C\ for cosmic microwave background anisotropics. 
IV. COUPLINGS TO NORMAL MATTER AND INSTABILITY 

Excitations of the scalar field from the model in the previous section have negative energy. 
The existence of negative energy particles may cause the system to be unstable due to 
interactions involving these particles. 

In the next section we examine this possibility by estimating the tree level decay rate 
of negative energy phantom particles into other phantoms and gravitons. Before delving 
into that calculation, we will first try to build some intuition about possible instabilities 
by considering the classical evolution of a simple system, a coupled pair of simple harmonic 
oscillators. Conservation of energy limits the phase space of a coupled pair of oscillators when 
both oscillators have positive energy. Though the oscillators may exchange energy, neither 
oscillator may ever reach an energy greater than the total initial energy of the system. 
Allowing one of the oscillators to have negative energy removes this limitation on the phase 
space. In this case the positive energy oscillator may increase its energy to any level so long 
as the negative energy oscillator decreases its energy by a compensating amount. Hence 
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the positive energy oscillator may reach arbitrarily high energies, and the negative energy 
oscillator may reach correspondingly large negative energies. 

While conservation of energy does not prevent the oscillators in this simple model from 
reaching arbitrarily large energies, there is no guarantee that the system will be unstable in 
all regimes. In fact, the following analysis will show that for a weak coupling, the evolution 
of the energy of each oscillator exhibits a stable oscillation for an arbitrarily long time. 

We consider a coupled pair of oscillators, one ip with positive energy, representing normal 
matter, and the other <f> with negative energy, representing the phantom component. 

1 



1 • 

2 



(22) 



where A is a dimensionless coupling constant. We are interested in the evolution of the 
energy of each oscillator, 

\i> 2 + m\ 



( \tf + ml 



The equations of motion are 



r = - 



1 - A^ 2 1 , 



where we have defined dimensionless rescaled variables 



(23) 



(24) 



A 



<j>/M 
i/j/M 



(25) 



In this section a prime denotes diferentiation with respect to the dimensionless time param- 
eter r, and M is the scale of the initial displacements of the oscillators. 

We explore the stability of this model by integrating the equations of motion (|24|) nu- 
merically for a range of parameters and initial conditions. Of particular interest are the 



16 



simulations with = 0, in which one may think of the ip oscillator as analogous to the 
massless graviton. Three such simulations are shown in Fig. |7| In each of these integrations 
the oscillators are started at rest, displaced a distance M from the origin. These plots are 
similar to those obtained for other initial conditions, and allowing the if) field to have a mass 
does not qualitatively change the plots shown here. Note that for A < 1, p^ and p^ exhibit 
a stable oscillatory behavior, while for A > 1, p^ and — p^ rapidly grow. From this analysis 
we conclude that this simple model is stable for small enough coupling but exhibits an in- 
stability for a large coupling. While the analysis of this simple model does not prove that a 
phantom field in the cosmological context will be stable to excitations and interactions with 
the gravitational field, it does provide some evidence that stability is possible. 

To make some connection to the cosmological model of the previous section, we consider 
how small the coupling in this simple model should be in order to have a stable solution. 
If we imagine a period of early universe inflation, then one expects perturbations in the 
phantom and gravitational fields to be of order M ~ 10 _5 Mp. The mass squared of the 
phantom field is 

<=%~wr (10_33eV)2 • (26) 

The restriction A < 1 then implies that the coupling in the Lagrangian satisfies A < 10~ 110 . 
While this at first seems to be an absurdly small number, we will see in the next section 
that, given the cosmological constraints of the previous section, such a coupling naturally 
results from considering perturbations of the phantom and gravitational fields. 



V. PHANTOM DECAY RATE IN FIELD THEORY 

The toy model of the previous section demonstrated that coupling a phantom oscillator 
to an ordinary oscillator results in a system which may or may not be unstable, depending 
on the magnitude of the coupling between them. It is far from clear, however, that this 
conclusion extends immediately to field theory. Roughly speaking, oscillators with some 
frequency correspond to field-theory modes of fixed wavelength; even if the model is stable 
when only certain wavelengths are considered, it does not follow that stability continues to 
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FIG. 7: Evolution of the energies of a coupled pair of oscillators, one with positive energy, 
(thin line), and one with negative energy, p^ (thick line), as a function of the rescaled coupling A. 
The oscillators are both started at rest displaced a distance M from the origin, and the energies 
are expressed in units of (m^M) 2 . Note that for A < 1, p^p and —p^ exhibit a stable oscillatory 
behavior, while for A > 1, p^ and — p^ rapidly grow. 

obtain when integrating over all momenta. 

One way of stating this concern is to ask about the decay rate of single particles that 
would conventionally be stable. Because excitations of the phantom field have negative 
energy, we could imagine a single particle decaying into a large number of phantoms and 



18 



ordinary particles. The rate for this process can be calculated (at tree level) using ordinary 
Feynman diagrams. We will find that the rate is infinite when arbitrarily high momenta 
are included, but can be rendered finite if a cutoff (which might arise, for example, from 
higher-derivative terms in the action) is introduced. 

Since we are interested in the allowed decay modes of phantom particles and their 
associated decay rates, it is instructive to first analyze the kinematics of reactions involv- 
ing phantom fields. Let us adopt the convention that ordinary particles have 4-momenta 
with positive timelike component. Then, since our field has a negative kinetic term, the 
corresponding 4- momentum will have a negative timelike component. This leads to the fol- 
lowing useful dictionary for translating between kinematically allowed reactions for phantom 
particles and those for ordinary particles. 

Denote ordinary particles as ip. To ask whether a certain reaction is kinematically allowed, 
we can just switch the phantom particles from the right side to the left and vice- versa, and 
ask whether the resulting reaction would normally be allowed. Consider for example the 
decay of an ordinary particle into another particle plus a phantom: 

V>i — > ^ 2 + • (27) 

This will be allowed if the reaction 

V>i + — > ^ 2 (28) 

would be allowed by conventional kinematics — in particular, if the mass of ordinary particle 
2 were greater than the sum of the masses of ordinary particle 1 and the phantom. So it is 
clear that ordinary particles can decay into heavier particles plus phantoms. For example, if 
the electron were coupled to the phantom field, processes such as 

e~ — > fi~ + v e + + (29) 

would be allowed. The muon would then decay back into an electron, as part of a potentially 
disasterous cascade. 

Next we turn to the decay of phantoms. First consider decays into ordinary particles 

— > Vi + ^2 (30) 
19 



will be allowed if 

|O)^0 + ^i + ^2 (31) 

would be conventionally allowed, which it is not. However, consider a decay into one phantom 
and one ordinary particle: 

01 — + 02 + V (32) 

will be allowed if 

02 — > 01 + ^ (33) 

would ordinarily be, which requires > + m^. So a phantom can decay into a heavier 
phantom plus a not-too-heavy ordinary particle. 

If there is only one kind of phantom, one might think it would be stable since there would 
be no heavier phantoms to decay into. However, several lighter particles can mimic the 
4-momentum of a heavier particle. Consider the decay of one phantom into two phantoms 
plus an ordinary particle: 

01 — > 02 + 03 + $ • (34) 

This will be allowed if 

02 + 03 ► 01 + i> (35) 

would ordinarily be, which it aways is for large enough relative velocities of 02 and 03. 

In summary, if there is only one kind of phantom particle (with a unique mass), it can 
only decay by emitting at least two more phantoms, plus at least one ordinary particle. 
Ordinary particles, meanwhile, may decay into phantoms plus other ordinary particles with 
a larger effective mass than the original. A special case to these rules comes from massless 
particles. Although it would seem kinematically possible for massless particles to decay into 
massive particles plus phantoms, massless particles cannot decay in flat space (a reasonable 
approximation in the backgrounds we are considering) simply because there is no rest frame 
in which to calculate the rate — no proper time elapses along a null path. 

Like any other dark-energy scalar field, the phantom should be weakly coupled to ordinary 
matter (or it would have been detected through fifth-force experiments or variation of the 



constants of nature 



33j). We therefore restrict our attention to only gravitons and phantoms. 
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FIG. 8: Decay of a phantom particle into two other phantoms and a graviton 

With the above rules in mind, the phantom decay channel involving the smallest number of 
particles is 

(pi -> h + 0i + 02 , (36) 

where h is a graviton and is a phantom particle, illustrated in figure |H1 

Consider the specific model investigated in section ITTT1 a phantom scalar with potential 

V(4>) = Vbe-W . (37) 

We will first consider this potential expanded as a power series around some background 
value 0o ~ Mp. Gravitons, meanwhile, may be represented by transverse-traceless metric 
perturbations g^ u = gffj + h^ u , where gffj is the background Robertson- Walker metric. 
Strictly speaking, there is no way for a single graviton to couple non-derivatively to the 
potential, simply because there is no way to construct a scalar from a single traceless h^. 
But to avoid a surfeit of indices, and because it won't affect the final answer, we will simply 
think of the graviton as a dimensionless scalar field h; to get a canonically-normalized field 
we multiply by Mp. 

To study the decay (JHHJ), we require the interaction part of the Lagrangian to first order 
in h and third order in 0, which is 

Cj = ]f L(Mp/ i )iv"'(0o)0 3 



= A cfr (M P /i)0 3 , 
21 



(38) 



where 

A cff = 20 O ^M „ 10 -«> . (39) 

r Mp Mp v ; 

The decay rate of a phantom particle through this channel is 

r = ^/pfepfepfeMW^-«- ■ (40) 

where the matrix element |.M| is just A e fr at tree level. To get an upper limit on the 
reaction rate, and hence a lower limit on the timescale, we assume approximate isotropy, so 
d 3 p ~ \p\ 2 dp. Because the relevant momenta are very large (and the masses very small), we 
can also approximate E ~ p. Putting all of this together we obtain 

r ~ — / \Ph\dph I f \p0 2 \dp0 2 5 {4) (yty. -pfa -P4, 2 -p h ) ■ (41) 

If we take the limits on the integrals to be oo, the decay rate is clearly infinite. Hence, 
the answer to our investigation into instability seems very clear: the theory is dramatically 
unstable, as individual particles rapidly decay into cascades of phantoms and gravitons. 

However, our philosophy has been to think of this model as an effective theory valid at 
low energies. The reason why the decay rate diverges is because the phase space is infinite, 
since the phantoms can have arbitrarily large negative energies. Therefore, this result relies 
on taking the calculation at face value up to infinite momentum transfer. Instead, we should 
only trust the phase-space integrals up to some cutoff where new physics might enter. For 
example, we could imagine a higher-derivative term of the form 

L ~ -(c>0) 4 , (42) 

which would eventually dominate over the negative-kinetic-energy term we have already 
introduced. 

We have not investigated closely the properties of phantom models with higher-derivative 
terms. Instead, let us crudely approximate the effect of a momentum cutoff at a scale A by 
only allowing the phase-space integrals to range up to that cutoff. Looking at (|41|). we can 
estimate the truncated decay rate as 

A 2 



Ke— ■ (43) 



rritf, 
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The timescale for decay is r = Despite any fundamental instability, a model will be 

phenomenologically viable if the lifetime is greater than the Hubble time Hq 1 ~ 10 60 Mp . 
Recalling that ~ 10 _60 Mp and A c ff ~ 10~ 120 , the lifetime in units of the Hubble time is 

tfor~10 120 (^) 2 . (44) 

In other words, the lifetime from this decay channel exceeds cosmological timescales so long 
as A < 10 60 Mp, which is certainly not a stringent constraint. 

However, the infinite phase space for this one decay is not the only infinity we have to 
deal with; a single phantom can decay into arbitrary numbers of gravitons and phantoms, 
so we must sum over all the channels. For each new final-state particle F gets multiplied 
by a factor of (A/Mp) 2 . The A 2 comes from the new momentum integral, and the Mp -2 is 
introduced into A 2 ff . So, if we denote the decay rate (j4*3^) by r , then if there are n additional 
particles in the final state the total decay rate is given by 

/ A \ 2n 

r » = U) r »- < 45 > 

Therefore, the total decay rate is 

oo 

r totaJ = £(n + l)r n , (46) 

n=0 

where the factor (n + 1) comes from the number of different ways the final state can be 
composed of gravitons and phantoms. This yields 

(47) 

Thus, the decay rate remains of order To so long as A is not larger than Mp. This still seems 
like a comfortable result. 

There is, however, one more possibility to be accounted for. Part of the reason we obtained 
reasonable decay rates from interactions originating in the potential V((f>) was because of 
the small value of Vo ~ (10~ 3 eV) 4 . This suppression is lost if we consider couplings of 
gravitons to derivatives of 0. Since we are claiming that our model is simply an effective field 



total 



A 

MZ 
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theory, we are obligated to include all possible non-renormalizable interactions, suppressed 
by appropriate powers of the cutoff scale. Consider for example the operator 

C = J^X^r^Wd^ , (48) 

where (3 is a dimensionless coupling. (We will use the same cutoff A for our non- 
renormalizable terms as in our momentum integrals; the Mp in the denominator comes 
from normalization of h^ v '.) This term will also contribute to the interaction shown in Fig- 
ure |H1 Following similar logic as above (including two powers of the momentum in M. from 
the derivatives), we get 



Using Hq ~ m,/, ~ 10 _60 Mp, the lifetime in units of the Hubble time is now 



r ~ ^—T^y • (49) 



^~^(^) 4 . (50) 

Therefore, if (3 is of order unity, to obtain cosmologically viable decay rates (H t > 1) we 
require the cutoff to be 

A < fCr 30 M P ~ l(T 3 eV . (51) 

This is a much smaller cutoff than was required when we considered coupling through the 
potential, since the small prefactor Vo is not around to help us. 

The above result is alleviated somewhat by imposing an approximate global symmetry 
on the theory, that the Lagrangian density be invariant under 0^0 + constant. Such a 
symmetry is quite reasonable, as it is the only known way to ensure both a nearly-flat po- 
tential and appropriately small couplings to ordinary matter [33]. In this case, the irrelevant 
operator of lowest dimension leading to phantom decay into gravitons is 

C = m|a7 KMp/n^<R0] 2 . (52) 



The decay rate is then 



with an associated lifetime 



7 2 A 6 



(53) 
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which, if 7 ~ 1, leads to the relaxed requirement A < 10~ 20 Mp ~ 100 MeV. 

Under our most optimistic assumptions (of an approximate global symmetry), we there- 
fore find that the momentum cutoff characterizing our effective theory must be less than 
100 MeV to guarantee that the instability timescale is greater than the age of the universe. 
We find this value to be uncomfortably low, although not necessarily impossible; dark-energy 
models are typically characterized by energy scales ~ 10~ 3 eV (for the vacuum energy to- 
day) and ~ 10~ 33 eV (for the mass of the scalar field), so perhaps the scalar field theory 
is only valid up to relatively low momenta. Alternatively, one might imagine searching for 
some mechanism which would suppress derivative couplings, leaving only the couplings of 
gravitons to the potential, which were consistent with a cutoff as high as the Planck scale. 



VI. CONCLUSIONS 



There is no doubt that the discovery of a new component of the energy density of the 
universe has profound implications for the relationship between particle physics and gravity. 
Whether this component be a pure cosmological constant (whose magnitude we have no idea 



how to understand), a dynamical component 
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(whose special interactions give rise to the tiny vacuum energy we observe) or 
an as yet unimagined source, its nature is an outstanding problem of fundamental physics. 
The data thus far, although pointing conclusively to the existence of this dark energy, do 
not allow us to distinguish between competing scenarios. In particular, much of the allowed 
parameter space lies in the region w < — 1 in which some cherished notion behind our present 
theories must be sacrificed. 

There are several ways in which one may achieve w < — 1, including purely negative kinetic 
terms, non-minimal kinetic terms, and scalar-tensor theories. In addition, one might imagine 
that effective superexponential expansion of the universe might be obtained by modifying the 
Friedmann equation. While this is at odds with a purely four-dimensional general relativistic 
description of the universe, such an effect might be obtained in the context of brane-world 
models [44]. However, any such modification of the Friedmann equation must avoid conflict 



with the precision predictions of primordial nucleosynthesis j45, 



46, 
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In this paper we have taken seriously the possibility that w may be less than —1 and 
have asked the question: "From a particle-physics and general-relativistic point of view, can 
such theories be made consistent?" In particular, we have considered a specific toy model in 
which the null dominant energy condition is violated and hence the resulting space-time may 
be unstable. In this model the cosmology is well-behaved and the theory may be constructed 
so that it is stable to small, linear perturbations. When we consider higher order effects, 
however, the model fails to remain stable. The central result of this paper is a field theory 
calculation of the decay rate of phantom particles into gravitons. This decay rate would 
be infinite if the phantom theory was fundamental, valid up to arbitrarily high momenta, 
and would render the theory useless as a dark energy candidate. We therefore consider the 
phantom theory to be an effective theory valid below a scale A. Interestingly, couplings of 
gravitons to an appropriate scalar potential do not lead to decay of phantoms into gravitons 
and other phantoms on sub-Hubble timescales, so long as the cutoff A is below the Planck 
scale. However, in such an effective field theory approach, we are mandated to include in the 
Lagrangian operators of all possible dimensions, suppressed by suitable powers of the cutoff 
scale. In particular, we must include couplings of gravitons to derivatives of the phantom 
field. We find that such operators, even though they may be of high order, can lead to 
unacceptably short lifetimes for phantom particles unless the cutoff scale is less than 100 
MeV, so new physics must appear in the phantom sector at scales lower than this. 

Our analysis demonstrates that a model with a well-behaved cosmological evolution and 
stability to linear perturbations may still exhibit instability due to higher order interactions. 
We may therefore ask another crucial question: "Should observers seeking to constrain 
cosmological parameters take seriously the possibility that w < —1?" Unfortunately the 
answer is somewhat ambiguous. On the one hand, we know of no easy way to construct 
a viable model of this sort; on the other, it is certainly conceivable that new physics in 
the dark-energy sector kicks in at low scales to render a phantom model stable, or that 
phantom behavior is mimicked by even more exotic mechanisms (such as modifications of 
the Friedmann equation). Therefore, whether or not observations constrain w to be greater 
than or equal to —1 is still an interesting question, although there is a substantial a priori 
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bias against the possibility. For theorists, our conclusion is more straightforward: the onus 
is squarely on would-be phantom model-builders to show how any specific proposal manages 
to avoid rapid vacuum decay. 
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